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(N . Abstract 

We have proposed that the cosmic ray spectrum "knee" , the steepening of the cos- 
r** ■ mic ray spectrum at energy E £ 10 15 - 5 eV, is due to ' 'new physics", namely new 

interactions at TeV cm energies which produce particles undetected by the exper- 
' imental apparatus. In this letter we examine specifically the possibility that this 

interaction is low scale gravity. We consider that the graviton propagates, besides 
the usual four dimensions, into an additional 5, compactified, large dimensions and 
we estimate the graviton production in pp collisions in the high energy approx- 
imation where graviton emission is factorized. We find that the cross section for 
graviton production rises as fast as (^/s/Mf) 2+s , where Mf is the fundamental 
scale of gravity in 4 + 5 dimensions, and that the distribution of radiating a fraction 
y of the initial particle's energy into gravitational energy (which goes undetected) 
behaves as Sy 6-1 . The missing energy leads to an underestimate of the true energy 
and generates a break in the inferred cosmic ray spectrum (the "knee"). By fitting 
the cosmic ray spectrum data we deduce that the favorite values for the parameters 
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of the theory are Mf ~ 8 TeV and 5 = 4. 
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Cosmic rays (CR) is a subject almost a century old; despite the amount of 
knowledge accumulated since their discovery, there is a number of outstanding 
issues concerning their origin, acceleration and composition. One of the most 
interesting features is the breadth of their energy spectrum which extends over 
11 orders of magnitude up to and beyond 10 20 eV (see [1] for a recent review). 

The CR spectrum can be described as a broken power law E -1 , with 7 obtaining 
three different values in three different energy regimes: In the range 10 9 — 10 15 ' 5 
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eV, 7 ~ 2.75. At energies greater than 10 15,5 , the "knee", the spectrum steep- 
ens to 7 ~ 3 out to an energy E ~ 10 18 eV (the "ankle"), at which point it 
becomes flatter again, 7 ~ 2.2 — 2.5 and extends to ~ 10 20 ' 5 eV, beyond which 
the flux is too small to give a meaningful flux even with the largest current 
detectors. 

The currently most popular models for cosmic ray acceleration (supernova 
shocks) can account for the CR spectra up to energies 10 1A Z eV only [Z 
is the corresponding nuclear charge). Measuring the composition and under- 
standing the origin of the cosmic ray spectrum above these energies remains a 
challenge both experimentally and theoretically. Of particular interest in this 
respect is the origin of the steepening of the CR spectrum at the "knee" : As we 
suggested in earlier work [2], this (in fact any) spectral steepening cannot be 
accounted for by considering the spectrum to be the sum of two independent 
CR components, presumably the output of two different kinds of CR sources 
(an explanation that would be valid for a flattening of the spectrum such as 
that observed at the "ankle"). Furthermore, we suggested in the above refer- 
ence that the simplest explanation of the CR spectrum break at the "knee" is 
that it is due to a novel channel in the physics of high energy pp-collisions, 
as proposed also by others in the past [3]. We also noted there that the cm 
energy corresponding to the "knee" is a few TeV, a scale tantalizingly close 
to the energy at which the emergence of "new physics" is anticipated on very 
general grounds. We advanced the idea that the new interaction, operative 
at and above these energies, results in particles that are not detected by the 
associated experimental devices, leading to an underestimate of the incident 
particle's energy for energies above that of the "knee". As a result, a cosmic 
ray spectrum which is a single power law in energy, will develop an increase 
in its slope (a "knee") at the energy at which the new interaction turns-on, 
with the spectrum reverting to its original slope when saturation occurs. 

The main motivation for introducing new physics comes from the need to 
provide a unified theory in which two disparate scales, i.e. the electroweak 
My/ ~ 100 GeV and the Planck scale Mp ~ 10 19 GeV can coexist (hierar- 
chy problem). Recently a novel approach has been proposed for resolving the 
hierarchy problem [4] . Specifically, it has been suggested that our four dimen- 
sional world is embedded in a higher dimensional space with D dimensions 
of which 5 dimensions are compactified with a relatively large radius. While 
the Standard Model (SM) fields live on the 4-dimensional world (brane), the 
gravitons can propagate freely in the higher dimensional space (bulk). The 
fundamental scale Mf of gravity in D dimensions is related to the observed 
4-dimensional Newton constant Gn by 
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where V$ is the volume of the extra space. For a torus configuration 



V s = (2nR) 5 



(2) 



with R the common radius of the large extra dimensions. For a given 5, a 
sufficiently large radius R can then reduce the fundamental scale of gravity 
Mf to TeV energies which are not too different from M\y, thereby resolving 
the hierarchy problem. 

The D— dimensional graviton, when reduced to 4— dimensions, gives rise to 
spin 2 particles (the massless graviton and its massive Kaluza-Klein excita- 
tions), spin 1 particles and spin zero particles. The vector particles do not 
couple to the energy momentum tensor and are ignored, while the scalar par- 
ticles couple to the trace of the energy momentum tensor, giving negligible 
contributions to the high energy regime we are interested in [5,6]. Thus we 
take into account only the production of spin 2 gravitons in the pp-collisions 
we consider. 

We study the pp — > pp collision at the cm system, where each particle car- 
ries energy y/s/2. At high energies the dominant contribution to graviton (h) 
production (pp — > pph), originates from collinear bremsstrahlung of gravitons 
from each external line. In this configuration, graviton emission is factorized 
as the probability that a proton (incoming or outgoing) loses a fraction y of 
its energy via graviton radiation. The Kaluza-Klein (KK) excitations of the 
graviton have the same couplings to ordinary fields as their massless zero- 
mode, i.e. they couple to the 4— momentum of the fermion they are attached. 
Within this reasonable approximation of factorized emission, we obtain that 
the cross-section for graviton production in pp collisions, with the graviton 
carrying energy e, is given by 



where o"o is the pp — > pp cross-section (it rises slowly with energy), m p is the 
proton mass and g(e) is the multiplicity factor counting the number of KK 
gravitons contributing to the process. The mass of each KK mode corresponds 
to the modulus of its momentum in the direction transverse to the brane 



with n = (ni, n 2 , ...,ns). At fixed e, all gravitons with m n < e contribute, their 
number being the volume of a 5— dimensional sphere with radius n max = Re. 
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Thus 
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Substituting expression (5) into eq. (3), using relation (1), absorbing constants 
and 5— dependent factors into a rescaled Mf, we find that the cross-section for 
a fractional energy loss y (y — 2e/y/s) is 

-J^(PP -»■ PPh) = cr F(s)f(y) (6) 



with 
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f(y) = sy 
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We observe that the cross-section for graviton emission rises fast, as a power 
of the energy (i.e. oc (y/s / 'M /) 2+5 ) and while it is unimportant at energies 
below the scale Mf, it is the dominant process at and just above Mf. The dis- 
tribution of the energy radiated into massless gravitons is represented by the 
well-known infrared 1/y behavior, but the large number of KK modes converts 
this distribution into y 5 ^ 1 . We emphasize again that within our approxima- 
tive scheme we pick-up the leading contributions originating from collinear, 
soft, gravitational radiation. For the purposes of the present investigation we 
consider this approximation as well justified. 

In the lab system a cosmic ray particle of high energy hits a nucleon at rest 
in the Earth's atmosphere. For a fractional energy loss y in the cm system, in 
the lab system the corresponding event of total energy E' will be registered 
at the detector at an energy E — E'(l — y/2) 2 . Therefore, for a cosmic ray 
intensity 1(E), gravitational radiation will induce an inferred intensity N(E) 
of the form 

l 

N(E) = J dE'J dyf(y)I(E')5(E-E'(l-y/2f)) (9) 
o 



The cosmic ray interactions will proceed either through the standard channels 
with probability Pq(E), or through the production of gravitons with proba- 
bility Ph(E). We deduce that the observed cosmic ray intensity I Q b(E) at an 
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Fig. 1. The cosmic ray spectrum f(E) multiplied by E 3 for E > 10 14 eV. Long 
and short dashed lines are respectively the incident galactic and extragalactic 
components. Solid line is the combined spectrum expected to be measured for 
Mf = 8 TeV, 5 = 4. +'s are the Tibet data, diamonds the Casa Blanca data, 
squares the Fly's Eye data and x 's the AGASA data. 

energy E will be the sum of these two processes, i.e. 



I ob (E) = P I(E) + P h N(E) 



(10) 



where 



P o{E)= ^^=^_ (11) 

° K J a + a h 1 + F 1 ; 

PhiE)= ^t^ = TTF (12) 

a + a h 1 + t 



We assume that the incident cosmic ray spectrum 1(E) is composed of two 
components, the galactic component (Ig{E)) and the extragalactic compo- 
nent (Ie{E)) each of which is assumed to be a power law of indices 7 = 
2.75 and 2.2 respectively, with the galactic component being dominant at 
energies E 10 18 ' 5 eV and the extragalactic one at higher energies as dis- 
cussed in [2] (i.e. I G (E) oc E~ 2 - 75 exp(—E/E ), E ~ 10 18 - 5 eV and I E (E) oc 
E~ 2 - 2 exp(—E I 'Ex), E x ~ 10 20 - 5 eV). According to our suggestion, the observed 
spectrum I b(E) results from a migration of high energy points of energy E' 
to lower energy E [see Eq. (9)]. Therefore, the inferred spectrum is a sensitive 
function of the input spectrum 1(E). Indeed the resulting spectrum, shown in 
fig. 1, depends both on the parameters of the theory (Mf, 5) and the parame- 
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ters of the input spectrum, notably the value of the galactic cut-off energy E . 
Rather than varying freely the parameters of the input spectrum, we have cho- 
sen the values adopted in our previous work [2] for the reasons explained there. 
The uncertainty in the experimental data, exemplified by the disagreement be- 
tween two experiments in each energy range (i.e. Tibet - CASA BLANCA for 
E = 10 14 - 5 - 10 16 - 5 eV , Fly's Eye - AGASA for E > 10 17 - 5 eV) presented in 
figure 1, does not allow a high significance determination of the parameters of 
the theory. We found however that the data favor the values Mf ~ 8 TeV and 
S = 4 for the theoretical parameters, with the values S — 3, 5 being of lower 
significance as it was determined by the numerical study. Our deduced values 
respect all available bounds on Mf and 5. The most severe bounds come from 
the observations of SN 1987a, which constrain the rate of energy loss due to 
new particles [7]. Detailed caclulations [8] provide the limit Mf ^ 1 TeV for 
5 = A. 

Cosmic accelerators accelerate particles (the cosmic rays) to energies much 
higher than those achieved at terrestrial laboratories. As such they present 
us with the opportunity of detecting the signatures of processes operative at 
energies not yet accessible at terrestrical accelerators. As mentioned in the 
beginning of the present note, the lowering of the fundamental Planck scale 
Mf to TeV energies, implies that the gravitational interaction becomes strong 
at these energies. Any scattering process at cm energies E ~ Mf should be 
accompanied by abundant graviton production. Since the energy of the cosmic 
rays can by far exceed this threshold domain, graviton emission should be pre- 
eminent in CR interactions. The essence of our proposal is that the signature of 
emission of these gravitons, which go undetected in the cosmic ray air shower 
arrays, is the observed "knee structure" in the CR spectrum. 

We would like to conclude with a few general remarks. How reliable is our 
calculation? It is obvious that our cross section for graviton emission rises 
very fast and at some energy unitarity will be violated. The KK formulation 
we employed is an effective theory one and at some energy M s (M s is the 
string scale with M s > Mf) we have to resort to the underlying theory, a 
string theory. With gravity dominant at TeV scales, gravitational radiation 
is only one of the possible manifestations. Exchange of virtual KK gravitons 
enhances the cross-sections [9] and this enhancement might help to resolve 
the paradox with the CR events above the GZK cut-off [10]. It has been 
suggested also that black holes may in fact be produced in pp collisions [11]. 
Their subsequent decay through the Hawking radiation should lead to events 
of high multiplicity and large sphericity at the cm. It is of interest that detailed 
analysis of the EAS data suggest an apparent, very sharp change in the CR 
composition to almost exculsively Fe just above the "knee" (see fig. 5 of [12]). 
This apparent change is qualitatively of the form expected by a sharp increase 
in the interaction cross section and an ensuing dispersion of the available 
energy to a large number of secondary particles. The process described just 
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above might then provide an account of this fact, though we do not believe 
that this interpretation can be considered at this point as unique. Finally, on 
the cosmic ray physics side, our proposal implies that the galactic CR sources 
produce power law spectra extending to the "ankle" rather than the "knee" (as 
thought on the basis of acceleration in SN shocks). Given that the latter do 
constitute CR sources to these energies our arguments imply the presence of 
two independent galactic CR components: one due to SN shocks, while clues 
of the origin of the second, being the observed anisotropy of CR at E ~ 10 18 
eV toward the galactic center [13]. It appears that cosmic rays rather than 
being the messengers, constitute the message itself. 
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